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Abstract 

In this paper, we study the problem of power allocation for streaming multiple variable-bit-rate 
(VBR) videos in the downlink of a cellular network. We consider a deterministic model for VBR video 
traffic and finite playout buffer at the mobile users. The objective is to derive the optimal downlink 
power allocation for the VBR video sessions, such that the video data can be delivered in a timely 
fashion without causing playout buffer overflow and underflow. The formulated problem is a nonUnear 
nonconvex optimization problem. We analyze the convexity conditions for the formulated problem and 
propose a two-step greedy approach to solve the problem. We also develop a distributed algorithm based 
on the dual decomposition technique. The performance of the proposed algorithms are validated with 
simulations using VBR video traces under realistic scenarios. 

Index Terms 

Convex optimization; distributed algorithm, downlink power control, video streaming, variable bit 
rate video. 

I. Introduction 

According to a recent study by Cisco, data traffic over wireless networks is expected to increase by a 
factor of 66 times by 2013. Much of the increase in future wireless data traffic will be video related, as 
driven by the compelling need for ubiquitous access to multimedia content for mobile users. Such drastic 
increase in video traffic will significantly stress the capacity of existing and future wireless networks. 
While new wireless network architectures and technologies are being developed to meet this "grand 
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challenge" lll-llll, it is also important to revisit existing wireless networks, to maximize their potential 
in carrying real-time video data. 

Quality of service guarantee in wireless networks is a challenging problem that has attracted tremendous 
efforts ll6l-|[lTl. In this paper, we consider the problem of streaming multiple videos in the downlink of a 
cellular network. The system is interference limited: the capacity of a specific mobile user depends on the 
Signal to Interference-plus-Noise Ratio (SINR) at the user, which is a function of the power allocation 
for all the mobile users. Therefore, effective downlink power control is necessary for such a wireless 
video system to minimize the intra-cell interference for concurrent video sessions. 

We consider the challenging problem of streaming concurrent variable-bit-rate (VBR) videos in the 
cellular network. This is motivated by the fact that VBR video offers stable and superior quality over 
constant bit rate (CBR) videos. Furthermore, many stored video content are VBR. It is important to 
support such stored VBR videos over existing wireless networks without the need for transcoding. The 
main challenge in supporting VBR video stems from its high rate variability and complex autocorrelation 
structure, making it hard for network control and may cause frequent playout buffer underflow or overflow. 
In this paper, we adopt a deterministic traffic model for stored VBR video, which jointly considers frame 
size, frame rate, and playout buffer size ||T2| - ||T41 . Unlike prior work, we exploit effective downlink 
power control to adjust the downlink capacities based on prior knowledge of frame sizes and palyout 
buffer occupancies. 

Specifically, we present a downlink power control framework for streaming multiple VBR videos in a 
cellular network. With the deterministic VBR video traffic model, we formulate an optimization problem 
that jointly considers donwlink power control, intra-cell interference, VBR video traffic characteristics, 
playout buffer underflow and overflow constraints, and base station (BS) peak power constraint. The 
objective is to maximize the total throughput, which can achieve high playout buffer utilization. As a 
result, playout buffer underflow or overflow events can be minimized. We analyze the convex/concave 
regions of the formulated problem and develop a two-step downlink power allocation algorithm for 
solving the problem. We also develop a distributed algorithm based on the dual decomposition technique 
from convex optimization, in order to reduce the control and computation overhead at the BS. We 
evaluate the performance of the proposed distributed algorithm with simulations using VBR video traces. 
Our simulation results verify the accuracy of the analysis and demonstrate the efficacy of the proposed 
algorithms. 

The remainder of this paper is organized as follows. The system model is presented in Section JI] 
We develop a two-step algorithm to solve the power allocation problem in Section Jill and a distributed 
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algorithm based on dual decomposition in Section JV] Simulation results are presented in Section fV] and 
related work are discussed in Section |VI] Section I VII I concludes this paper 



II. Network and Video System Model 

We consider the downlink of a cellular network. In the cell, a BS streams multiple VBR videos 
simultaneously to mobile users in the cell, which share the downlink bandwidth. We assume the last-hop 
wireless link is the bottleneck, while the wired segment of a session path is reliable with sufficient 
bandwidth. Thus the corresponding video data is always available at the BS before the scheduled 
transmission time. 

VBR video traffic exhibits both strong asymptotic self-similarity and short-range correlation. A stochas- 
tic model capturing the complex auto-correlation structure often requires a large number of parameters, and 
is thus hard to be incorporated for scheduling real-time video data. To this end, we adopt a deterministic 
model that considers frame sizes and playout buffers Let Dn{t) denote the cumulative consumption 
curve of the n-th mobile user, representing the total amount of bits consumed by the decoder at time t. 
The cumulative consumption curve is determined by video characteristics such as frame sizes and frame 
rates. Assume mobile user n has a playout buffer of size 6„ bits and its video has T„ frames. We can 
derive a cumulative overflow curve for the user as 

Bn{t) = m.\n{Dn{t - 1) + 6„, I?„(T„)}, < t < T^. (1) 

Bn{t) is the the maximum number of bits that can be received at time t without overflowing user n's 
playout buffer. Finally we define cumulative transmission curve Xn{t) as the cumulative amount of 
transmitted bits to user n at time t. To simplify notation, we assume the video sessions have identical 
frame rate and the frame intervals are synchronized. A time slot t is equal to the t-th frame interval, 
denoted as r. 

The three curves for user n are illustrated in Fig. [T] A feasible transmission schedule will produce a 
cumulative transmission curve Xn{t) that lies within Dn{t) and Bn{t), i.e., causing neither underflow 
nor overflow at the playout buffer. In practice, Dn{tys are known for stored videos and can be delivered 
to the BS during session setup phase, and i3„(t)'s can be derived as in ([T]). 

We consider subscribers in the cell and let U denote the set of users. In each time slot t, the BS 
transmits to each user n with power Pn{t) and the power allocation is P{t) = [Pi{t), ■ ■ ■ , Pn{t)]^. We 
also consider a maximum transmit power constraint P, i.e., J2neu Pn{t) < P, for all t. When the power 
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allocation P{t) is detemiined, the SINR at user n can be written as lITSl . |[T6l 

where P„ is the power allocated to user n, G„ is the path gain between the BS and user n, rjn is the 
noise power at user n, Ln is a constant for user n (e.g., processing gain), and /3 denotes the orthogonaUty 
factor, with < /3 < 1. In this paper, we consider the case /3 = 1, where the SINR of a user not only 
depends on its own power allocation but also the power allocations of other users. 

We assume slow-fading channels such that the path gains do not change within each time slot [15]. The 
downlink capacity C„(t) depends on the SINR at user n, the channel bandwidth B^, and the transceiver 
design, such as modulation and channel coding. Without loss of generality, we use the upper bound as 
predicted by Shannon's Theorem: 

Cn{P{t)) = log (l + ^n{P{t))) . (3) 

In time slot t, Cn{t)T bits of video data will be delivered to user n. The cumulative transmission curve 

Xn{t) is 

X„(0) = 0; Xn{t) = Xn{t - I) + Cn{t)r. (4) 
For a feasible power allocation, the cumulative transmission curves should satisfy 

Dn{t) < Xn{t) < Bn{t), for all n,t, (5) 
i.e., without causing playout buffer underflow or overflow. 
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From (|3]l~(l5]l, the lower and upper limit on the feasible SINR at user n can be derived as 

= max{exp{ -i°'^-ff;^-(^-^)i },7f} 

where 7*^* is the minimum SINR requirement imposed by the transceiver design. 7^*"(t) is the SINR 
that the just empties the buffer at the end of time slot t, without causing underflow; 7^'^^(t) is the SINR 
that just fills up the buffer at the end of time slot t, without causing overflow. 

Generally, feasible power allocation P{t) is not unique for a given set of VBR video sessions. Among 
the set of feasible solutions, a schedule that transmits more data is more desirable since it provides more 
flexibility for optimizing future power allocations. We formulate the problem of optimal downlink power 
control for VBR videos, termed Problem A, as 

(A) maximize EnGW log(l + 7„(t)) (7) 
subject to: 

ln{t) = ^ p u^, — ' for all n (8) 

(t ) < 7n (t ) < 7^"" it) , for all n (9) 

EneuPn < P. (10) 

In Problem A, the objective is to achieve the maximum buffer uitilization at the users, under playout 
buffer underflow and overflow constraints and BS maximum transmit power constraints. This is a nonlinear 
nonconvex problem, to which traditional convex optimization techniques cannot directly apply. Due to 
the large variability of VBR traffic, the SINRs may assume values ranging from very low to very high, 
to avoid playout buffer underflow and overflow. Thus the existing high SINR approximation HI 71 and 
low SINR approximation fTSi techniques cannot be directly applied. 



III. Two-Step Downlink Power Allocation 

In Problem A, we consider an interference-limited system, where the capacity of downlink n depends 
on the power allocations for all the users. In the following, we first derive conditions for the optimal 
solution, and then present a two-step power allocation algorithm for solving Problem A. 

Lemma 1: If there exists a feasible power allocation P{t) that achieves 7^"^(t) for all n, the solution 
is optimal. 
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Proof: If the feasible power allocation P{t) achieves 7™"^(t) for all n, then all the user buffers 
are full at the end of the time slot, according to (|6]l. The objective value ^ cannot be further improved 
without causing buffer overflow. Thus the solution is optimal. ■ 
Lemma 2: If the upper limit 7™"^(t) cannot be achieved for every user n, then the optimal power 
allocation P{t) satisfies J2neu Pn{t) = P- 

Proof: Consider a feasible power allocation P'{t) = [P{{t) , P^{t) , ■■■ , P'j^{t)]'^ and Y.n&u Pni^) < 
P. We can construct another feasible power allocation P"{t) = [Pi'{t),P^'{t),- ■ ■ ,P^{t)f, such that 
= K ■ P^{t), for all n, and k • YnauKit) = UnrnKi't) < P' where k > 1. For the SINR at 
user n, we have 

7„(P"(t)) = LnGnPnjt) 

J2kj^n GnPk{t) + Vn 
_ KLnGnP'iiit) 

Y.k^n K-GnP'k{t) + r]n 

KLnGnP'nit) 
Hk^n I^GnPkit) + Kljn 

= ln{P'{t)). 

It follows that Y.n£U log(l + ln{P"{t))) > En6W log(l + ln{P'{t))), sincc Iog(l + x) is an increasing 
function of x. 

Choosing k = P/ J2neu Pni'^)^ we can construct a feasible solution P"'{t) = k ■ P'{t), such that 
EneuPnit) = P.Thenwehave7„(P"'(t)) > jn{P'{t)) and EnG^/ log(l+7n(^^"'(t))) > E„ewlog(l+ 
jn{P'{t)))- That is, any feasible solution with J2neU P'ni^) < P will be dominated by feasible solutions 
with J2n€U Pni't) = P- We conclude that the optimal solution P{t) must satisfy J2neu Pn{t) = P- ■ 

We have the following result for the optimal solution of Problem A, which directly follows Lemmas [T] 
andlH 

Theorem 1: A solution to Problem A is optimal if (i) it achieves the maximum SINR 7^°^(t) for all 
n; or (ii) its total transmit power is P. 

Theorem [U implies that we can examine the SINR (or buffer) constraints and the peak power constraint 
separately. In the rest of this section, we present a two-step power allocation algorithm for solving Problem 
A. We first examine Problem A under condition (i) in Theorem [T] to obtain Problem B as 

(B) 7r^(t) = ^"^"f"|f , for all n, (11) 

Efc^n GnPk{t)+Vn 

subject to: 

EneuPn < P- (12) 



TABLE I 

Two-Step Power Allocation Algorithm: Step I 



1 BS obtains D„, and Bn, and computes 7™"^ for all user n; 

2 BS tests the existence of feasible solutions using (Tsj; 

3 IF J13b is solvable, compute its solution P{t); 

ELSE, go to Step II of the algorithm, as given in Table HH 

4 IF J2„eu ^"(*) ^ stop with the optimal solution P{t); 
ELSE go to Step II of the algorithm, as given in Table HH 



In Problem B, (fTTI) is a system of linear equations of power allocation P{t). Rearranging the terms, 
we can rewrite ([TTI l in the matrix form as: 

(I - F) P{t) = u, for P{t) y 0, (13) 

where I is the identity matrix, F is a x matrix with 

0, \i n = ra 

(14) 

^max jj^^^ Otherwise, 

and u = [r?i7rVL„Gi, r/27r"/inG2, • • • , VniTVLuGnV. 

Since all the variables are nonnegative, F is a non-negative matrix. According to the Perron-Frobenius 
Theorem, we have the following equivalent statements |[T9l : 

Fact 1: The following statements are equivalent: (i) there exits a feasible power allocation satisfying 
([T3] ): (ii) the spectrum radius of F is less than 1; (iii) the reciprocal matrix (I — F)~^ = J2T=o (F)*^ 
exists and is component-wise positive. 

Based on Theorem [T] and Fact[T] we derive the first step of the two-step power allocation algorithm, as 
given in Table Jl If Problem B is solvable, the Step I algorithm in Table U produces the optimal solution 
for Problem A according to Theorem [1] Otherwise, we derive Problem C by applying Lemma [2j as 

(C) maximize E„gw log(l + 7„(t)) (15) 
subject to: 

< P„(t) < for all n (17) 

EneuPnit) = P, (18) 
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where An = rjn/Gn is the ratio of noise power and channel gain, representing the quality of the user n 
downlink channel. and P^°'^{t) are solved from dD and (IT6] l. as 

p--(i) = ^--(P + A„)/(L„ + 7--) ^^^^ 

Since the total transmit power is P, the objective value in ([TSll and the SINR in (fT6l ) for each user 
only depends on its own power Note that all the constraints are now linear To solve Problem C, we 
examine the objective function to see if it is convex. We omit time index t in the following for brevity. 

Lemma 3: The capacity of each user n, Cn, has one inflection point P*: when P„ < P*, C„ is in 
concave; when Pn > P*, Cn is convex. 

Proof: Taking the first and second derivatives of the objective function ( fTSl ) with respect to P„, we 

have 

dCniPn) Ln{P + An) 



(20) 
(21) 

Since Pn < P and An > 0, both the first and second derivatives exist. Letting — gp 2" = 0, we derive 



dPn {P - Pn + An)[P + {Ln - l)Pn + An] 

d^CniPn) -Ln[{Ln " 2)(P + An) + 2(1 - + An) 



dPn^ [{P-Pn + An)'' + LnPn{P-Pn + An)Y 

An>0 

the unique inflection point 



When Pn< P*, it can be shown that ^'^^^f"^ < 0; when Pn > P*, it can be shown that ^^gp^^"^ > 0. 

■ 

The normalized capacities for a two-user system is plotted in Fig. |2l with the inflection points marked. 
It can been observed that the curves are concave on the left hand side of the inflection points and convex 
on the right hand side of the inflection points. The processing gain is usually large for practical systems 
(e.g., Ln = 128 in IS-95 CDMA). We assume L„ » 1 in the following analysis. 

Theorem 2: For Problem C, there can be at most two links operating in the convex region if Ln > 

(4P + 6An)/{P + 3An). 

Proof: The reflection point is P* = j^ff^^{P+An). As L„ oo, we have P* = 0.5(P+A„). Only 
one link can operate in the convex region due to constraint (1181 ). Since > 0, P* is an increasing 
function of Ln- When 1 < L„ < cx), we have P* < 0.5 ■ {P + An). Letting SP* = P, we have 

Ln = (4P + QAn)/{P + 3A„). ■ 

For a clean channel where An ~ 0, L„ > 4 will guarantee at most two links operating in the convex 
region. The following results are on the impact of channel quality An = rjn/Gn. 




Theorem 3: For a given L„, the inflection point P* is an increasing function of An- For two links i and 
j with the same transmit power P, if Ai < Aj, we have Ci{P,Ai) > Cj{P,Aj) and ^'^'gp.'^''* |Pi=p > 

Qp^ \P^=P > 0. 

Proof: The first part can be easily shown by the first derivative of P* with respect to A^, which 
is = 2(L^-i) ^' ^ ^- "'"^^ second part can be easily shown by evaluating ([TS] ). ([T6] l. and 

(Unil. ■ 

Theorem |3] shows that, for two Unks in the convex region with the same initial power P, allocating 
more power to the link with better quality can achieve larger objective value than alternative ways of 
splitting the power between the two links (i.e., achieving the multi-user diversity gain). Based on the 
above analysis, we develop the second step of the power allocation algorithm for solving Problem C, as 
given in Table |lll In Table |II1 Lines 3 ~ 4 tests the feasibiUty of the power allocation. If the sum of the 
total minimum required power is larger than the BS peak power, there is no feasible power allocation 
and there will be buffer underflow. In this case, we select users with "good" channels for transmission 
and suspend the users with "bad" channels. 

The Step II algorithm checks the three possible solution scenarios for Problem C depending on the 
network status and video parameters: 

• All Unks operate in the convex region; 

• One link operates in the convex region and the remaining links operate in the concave region 
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TABLE 11 

Two-Step Power Allocation Algorithm: Step II 



Initialization : 

1 BS obtains 6„, _D„, and Bn for all user n; 

2 BS computes 7;""^, 7;""% and P*, for all n; 

3 BS computes the minimum required sum power Pmin = X/new 
and gap Ap = P - Pmin; 

4 IF Pmin > P, remove links from W, according to descending order 

of An, until Pmin < P', 

5 Compute = for all P^^ > P- 
Phase 1 : 

6 Select all the users satisfying P""' < P* as a set U' C W; 

7 Solve Problem C under constraints 

Pr" < Pn < min (P^"^P:) and P^ < P' = P ^ EnsM' ■P-"'' 

where ZY' is the complementary set of U' , and obtain solution Pi; 

8 Calculate P„ by updating P™*" to the solution in Line 7 and assign 
the remaining power to the nodes in set U, in descending order of Rn\ 

9 Obtain the Phase 1 solution, Pp^ , and objective value fp^ ; 
Phase 2 : 

10 Select the link with the maximum R„, and assign all the available 
power P — Pmin to the link, until either all the power is assigned or 
the link attains power P™"""; 

11 IF there is still power to allocate, THEN select all the nodes in set U\n 
and repeat Lines 5 ~ 8; 

12 Obtain the Phase 2 solution, Pp^, and objective value fp^; 
Phase 3 : 

13 Select the first 2 links with the largest Pn's, and assign all the availble 
power P — Pmin to the links, until all the power is assigned or the links 
attains power P™""^, and repeat Line 11; 

14 Obtain the Phase 3 solution, Ppg , and objective value fp^ ; 
Decision : 

15 Choose the largest objective value among /p^, fp., and /p,, and stop 
with the corresponding power assignment; 



• Two links operate in the convex region and the remaining links operate in the concave region. 
Each of the three phases in Table JI] considers the optimality condition for one of the three scenarios. 
In particular, Phase 1 first optimizes the power allocation in the concave region and then allocates the 
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remaining power to the links that could be moved to the convex region. Phase 2 allocates as much power 
as possible to the link with the best quality, which could work in the convex region. Phase 3 attempts 
to move the second best link to the convex region if the total power constraint is not violated. Usually 
when Ln and n are large, Phase 3 will rarely occur due to the peak power constraint. 

In Table Ull Line 7 presents a convex optimization component, for which several effective solution 
techniques can be applied. In the following section, we describe a distributed algorithm for Line 7 based 
on dual decomposition. 



IV. Distributed Algorithm 

As discussed in Section Jill the core of the Step II algorithm is to solve Problem C in the concave 
region (see Fig. |2]|. In this section, we present a distributed algorithm for this purpose, where the users are 
involved in power allocation to reduce the control and computation overhead on the BS. In the concave 
region, we have Problem D as 



(D) maximize 
subject to: 



EneU log(l +7n(t)) 
LnP-nit) 



- , for all n 



< Pn{t) < min{P™'^^,P*},for all n 



EneuPnit) < Pt 



tot ) 



(23) 

(24) 
(25) 
(26) 



where Ptot < -P is the total power budget for the links in the concave region. For brevity, we define 
P^^ = minjP^'^^', P*} and drop the time slot index t in the following analysis. 

Introducing non-negative Lagrange multiphers A„, fin, and u for constraints (l25T l and (|26] |. respectively, 
we obtain the Lagrange function as 



C{P,X,fl, v) 



(27) 



log 1 + 



I-'nPn 



P Pfi + 



+ An(Pn — P^ 



+ 



EneW f^niPn^ - Pn) +V {Ptot - Y^neuPn) 

J2n£U ^n{Pn,Xn,tJ-n,l^) + {tJ-nPn^-XnP^^"') +J^Ptot, 



where 



jO-n{Pn, K,fJ'n, J^) = log 1 + 



Pn.Pn 



P ~ Pn + 



+ (An - /in - l')Pn- 



(28) 
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Since Cn only depends on user n's own parameters, we have the dual decomposition for each user n. 
For given Lagrange multipliers (or, prices) An, /in, and z>, we have the following subproblem for each 
user n. 

Pn{Kilin,i') = argmax £n(-Pn, A„, /U„, i>), for all n. (29) 

Subproblem (|29l ) has a unique optimal solution due to the strict concavity of £„. We use the gradient 
method |[20| to solve ( |29b . where user n iteratively updates its power Pn as: 

Pn{l + 1) (30) 
= [Pn{l) + 9{l)VnCn{Pn)\* 

where [•]* denotes the projection onto the range of The update stepsize 6 {I) varies in 

each step / and is determine by the Armijo Rule [20 1. Due to the strict concavity of the series 
{P„(l), Pn(2), • • •} will converge to the optimal solution Pn as / ^ oo. 

For a given optimal solution for problem (|29) . P = [A, • • • ,PnY, the master dual problem is as 
follows: 

minimize C{P,\,fl,v) (31) 

subject to: A„, > 0, for all n. (32) 



Since the objective function (yTjl is differentiable, we also apply the gradient method to solve the master 
dual problem EOl . where the Lagrange multipliers are iteratively updated as 

An(/ + 1) = [An(/) - ax{l) ■ for all n 

/xn(/ + 1) = [/.n(0 - ■ ^^§f^]+, for all n (33) 
_ z.(/ + l) = KO-a.(0-^^|^]+, 
where [•]+ denotes the projection onto the nonnegative axis. The update stepsizes are also determined by 
the Armijo Rule ll20l . As the dual variables A(/), /2(Z), z^(Z) converge to their stable values as / — )• cxd, the 
primal variables P will also converge to the optimal solution ||2TI . 

The distributed algorithm is given in Table |IlIJ where the above procedures are repeated iteratively. 
The BS first broadcasts Lagrange multipliers to the users. Each user updates its requested power as in 
(l30l) . using local information P^*", P^°-^ , P*, An, Ln, and BS peak power P. Each user then sends its 
requested power back to the BS, and the BS will updates the Lagrange multipliers as in (|33T l. And so 
forth, until the optimal solution is obtained. 
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TABLE III 

Distributed Power Control Algorithm 



1 BS sets / — and prices A„(?), 1/(1) equal to some 
nonnegative initial values for all n; 

2 BS broadcasts the prices to the selected users; 

3 Each user locally solves problem l |29t as in ( 1301 to obtain its 
requested power; 

4 Each user sends its requested power to the BS; 

5 BS updates prices A„(Z), /i„(i), as in i33\ and broadcasts 
new prices A„(i + + 1), u{l + 1) for all n; 

6 Set 1 = 1 + 1 and go to Step 3, until the solution converges; 



V. Simulation Results 

We evaluate the proposed algorithms with simulations, using a cellular network with 20 users. The 
downlink bandwidth is 1 MHz. The path gain averages are G„ = d~^, where dn is the physical distance 
from the BS to user n. The downlink channel is modeled as log-normal fading with zero mean and 
variance 8 dB IfTSII . The processing gains are set to L„ = 128 for all n. The distance dn is uniformly 
distributed in [100m, 1000m]. The device temperature is Tq = 290 Kelvin and the equivalent noise 
bandwidth is = IMHz. The BS peak power constraints is set to P = 10 Watts. We use three VBR 
movies traces. Star Wars, NBC News, and Tokyo Olympics, from the Video Trace Lib Il22l . Each playout 
buffer is set to 1.5 times of the largest frame size in the requested VBR video. 

In the simulations, the proposed power allocation algorithm is executed at the beginning of each time 
slot. In Fig. [3j we plot the cumulative consumption, overflow and transmission curves for NBC News 
transmitted to user 2. The top sub-figure is the overview of 10, 000 frames. We also plot the curves 
from frame 2, 620 to 2, 640 in the bottom sub-figure. We observe that the cumulative transmission curve 
X{t) is very close to the cumulative overflow curve B{t), indicating that the algorithm always aim to 
maximize the transmission rate as allowed by the buffer and power constraints. The playout buffers are 
almost fully utilized most of the time. There is no playout buffer overflow and underflow for the entire 
range of 10, 000 frames. Among the NBC News frames, frame 2, 625 is the largest frame. We let seven 
out of the 20 links playout this largest frame simultaneously at time slot 2, 625 in the simulation. There 
is no buffer underflow under such heavy load. 

In Fig. m we plot the power allocation and price updates for all the 20 links in one of the 10,000 
time slots. The power and prices converges in around 70 steps. The converged power vector is P = 
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Fig. 3. Transmission schedule for video NBC News to user 2. 
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Fig. 4. Convergence of power allocation and Lagrange multipliers. 



[0.0022, 1.396, 0.0356, 0.0024, 1.396, 0.0351, 0.0016, 1.396, 0.0356, 0.0026, 1.396, 0.0356, 0.0023, 
1.396, 0.0356, 0.0018, 1.396, 0.0356, 0.0034, 1.394] Watts. 

Finally, we compare the proposed algorithm with a diversity -aware power allocation scheme, where the 
BS allocates power according to channel quality. With this scheme, the best channel n will be assigned 
power to achieve its maximum required power Then the second best channel will be allocated 

power until its maximum required power is achieved, and so forth until all of P is allocated. 
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Fig. 5. Average playout buffer utilization. 



We simulate 50 users with the same network and video settings. We compare the algorithms by their 
average playout buffer utilization. In Fig. |5l we plot the average buffer utilization from frame 2, 000 to 
2, 999. It can be seen that the proposed algorithm consistently achieves high buffer utilization, ranging 
from 60% to 100%. The diversity scheme achieves buffer utilization around 30% except for frames from 
2,250 to 2,400. Such considerably higher buffer utilization translates to better video quality: there is 
no buffer overflow or underflow for proposed algorithm, while there is buffer underflow in 17% of the 
playout frames for the diversity scheme. 

VI. Related Work 

There have been several papers on VBR video over wired network. Due to long-range-dependent (LRD) 
VBR video traffic, the piecewise-constant-rate transmission and transport (PCRTT) method was used to 
optimize certain objectives while preserving continuous video playout. In lfT2l . Liew and Chan developed 
bandwidth allocation schemes for multiple VBR videos to share a CBR channel. In 1 13,1 , Salehi et al. 
presented an optimal algorithm for smoothing VBR video over a CBR link. Feng and Liu |[T4l introduced 
a critical bandwidth allocation algorithm to reduce the number of bandwidth variations and to maximize 
receiver buffer utilization. Due to the fundamental difference between wireless and wired links, these 
techniques cannot be directly applied to the problem of VBR video over wireless networks. 

The downlink power allocation problem was studied in lITSl . ||T6l . aiming to obtain the power allocation 
that maximizes a properly defined system utility. A distributed algorithm based on dynamic pricing and 
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partial cooperation was proposed. Deng, Webera, and Ahrens ll23l studied the achievable maximum sum 
rate of multi-user interference channels. These papers provide the theoretical foundation and effective 
algorithms for utility maximization of downUnk traffic, but the techniques used cannot be directly applied 
for VBR video over wireless networks with buffer and delay constraints. 

In ll24l . |[25i . the authors studied the problem of one VBR stream over a given time-varying wireless 
channel. In |[24l . it was shown that the separation between a delay jitter buffer and a decoder buffer is in 
general suboptimal, and several critical system parameters were derived. In ll25l . the authors studied the 
frequency of jitters under both network and video system constraint and provided a framework for quan- 
tifying the trade-offs among several system parameters. In this paper, we jointly consider power control 
in wireless networks, playout buffers, and video frame information, and address the more challenging 
problem of streaming multiple VBR videos, and present a cross-layer optimization approach that does 
not depend on any specific channel or video traffic models. 

VII. Conclusion 

We developed a downlink power allocation model for streaming multiple VBR videos in a cellular 
network. The model considers interactions among downlink power control, channel interference, playout 
buffers, and VBR video traffic characteristics. The formulated problem aims at maximizing the total 
transmission rate under both peak power and playout buffer overflow/underflow constraints. We presented 
a two-step approach for solving the problem and a distributed algorithm based on the dual decomposition 
technique. Our simulation studies validated the efficacy of the proposed algorithms. 
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